Magneto-optical rotation of spectrally impure fields and its nonlinear dependence on 

optical density 
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We calculate magneto-optical rotation of spectrally impure fields in an optical thick cold atomic 
medium. We show that the spectral impurity leads to nonlinear dependence of the rotation angle on 
optical density. Using our calculations, we provide a quantitative analysis of the recent experimental 
results of G. Labeyrie et al. [Phys. Rev. A 64, 033402 (2001)] using cold Rb 85 atoms.. 
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I. INTRODUCTION 

A very useful way to get important spectroscopic infor- 
mation is by measuring magneto-optical rotation (MOR) 
of a plane polarized light propagating through a medium 
0. Clearly it is desirable to obtain as large an angle of 
rotation as possible 0. 0] It is known that the angle 
of rotation is proportional to the density of the medium. 
Thus an increase in density will help in achieving large 
rotation angles. Recently, very large rotation angles in a 
cold sample have been reported Q. In this experiment, 
optical densities of the order of 10 — 20 were achieved. 
This experiment also reported a very interesting result, 
viz., the departure from the linear dependence of the ro- 
tation angle on the optical density. This departure has 
been ascribed to the nonmonochromatic nature of the in- 
put laser. The findings of the experiment warrant a quan- 
titative analysis of the dependence of the rotation angle 
on the spectral profile of the input laser. We present a 
first principle calculation of this dependence. Note that 
a nonlinear dependence on optical density cannot result 
from a simple argument based on the standard formula 
for the rotation angle 9: 



= nkl Re(x_ — x+) 



(1) 



where, k is the wave-number of the electric field, I is 
the length of the medium, and x± represent the linear 
susceptibilities of the medium for right or left circularly 
polarized components of the input field. Since x± are 
proportional to the number density, the rotation angle 
becomes proportional to the optical density a defined by 
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where, A = 27r/fc is the wavelength of the input field, N 
is the number density of the atomic medium. If one were 
to argue that spectrally impure laser field would replace 
X± by their averages over the width of the laser, then 8 
would continue to be proportional to the optical density 
a. A more quantitative analysis of the rotation angle is 
thus warranted. 

The organization of the paper is as follows. In Sec. ITT1 
we recapitulate the relevant equations for studies of MOR 



in atomic medium. In Sec. lIIII we describe a simple three- 
level atomic configuration in the context of MOR and 
discuss the effect of laser line-shape and magnetic field 
in a thick atomic medium. In Sec. IIVI we consider ex- 
perimental configuration used in £|. We show, how in 
an optically thick medium, one deviates from the linear 
dependence of rotation on optical density as a result of 
the spectral impurity of the laser field. We give a quan- 
titative analysis of the experimental data. 



II. BASIC EQUATIONS 

Let us consider that an atomic medium of length I is 
resonantly excited by a monochromatic i-polarized elec- 
tric field 



E(z, t) = x£e 



ikz — iujt 



+ C.C. 



(3) 



where, £ is the field amplitude and k — lu/c = 2n/\ is 
the wave number of the field, u> and A being the corre- 
sponding angular frequency and wavelength. The field is 
propagating in z-direction. Clearly, we can resolve the 
amplitude of the electric field into its two circular com- 
ponents as 



x£ = e + £ + + e_£_ 



(4) 



where, £± = £/v2 are the amplitude components along 
two circular polarization e± = (x ± iy)/y/2. 

While passing through the atomic medium these 
two circular components behave differently due to the 
anisotropy of the medium. Let \± be the susceptibilities 
of the medium corresponding to the two circular compo- 
nents. The electric field at the exit face of the medium 
can be written as 



ikz — iujt 
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where, we have assumed that the medium is dilute so 
that |47rx±| <C 1. In MOR, the polarization direction of 
the input electric field gets rotated due to the difference 
in their dispersions (phase shifts) in a non-attenuating 
medium. The electric field however, remains linearly po- 
larized after passing through the medium. In the present 
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case, because the atom interacts with near-resonant elec- 
tric field, the two circular components suffer attenuation 
(given by the imaginary part of x±) while propagating 
through the medium. Thus the medium, concerned here, 
is responsible for both dispersion and attenuation. We 
term the medium to be both circularly birefringent and 
circularly dichroic. The output electric field becomes el- 
liptically polarized under the action of such a medium. 
Thus to fully characterize the polarization state of the 
output field, one has to use the Stokes parameters 
The four Stokes parameters for an electric field are desig- 
nated by S a (a — 0, 1, 2, 3) and can be defined as follows: 
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where, I„ is the measured intensity along the polariza- 
tion direction n. Then the output polarization state can 
be characterized by the following three quantities: 



P 

tan 29 
tan 2(f) 



si 



si 



So 

lr (0 < < tt) , 
S 3 



S P 



(-7T/4 < 4> < 7T/4) , 



(7a) 
(7b) 
(7c) 



where, P is the degree of polarization, i.e., the ratio of 
the intensities of the polarized component to the unpo- 
larized one, 9 is the Faraday rotation angle of the input 
field and is measured between the major axis of the ellipse 
and the x-axis, <j) provides the ellipticity of polarization 
through the relation e = tan cf>. 

From Eq. JSJ one can express the output intensities 
along different polarization directions in the following 
way: 
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FIG. 1: Level diagram for a three-level configuration. The 
excited levels |±) (m e = ±1) are Zeeman shifted from the 
level m e = by an amount s. The detuning 5 is defined 
between the levels m e = and m g = 0. 



to be obtained from 



1 f°° 
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(9) 



For simplicity, we can adopt, say, a Lorentzian line-shape 
for the input 
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where, ui is the central frequency of the laser field and 
27 c is the full width at half maximum. Wc will demon- 
strate how the fluctuations of the input field leads to 
the nonlinear dependence of the rotation angle on opti- 
cal density. 



III. A SIMPLIFIED ATOMIC MODEL 

We first consider a three-level atom in V configura- 
tion (see Fig. ^) in order to uncover the effect of optical 
density on MOR. The levels |±) (J e = l,m e = ±1) are 
coupled to the ground state \g) (J g — 0, m g = 0) by 
two circular components a± of i-polarized electric field 
[Eq. J3J]. The excited level degeneracy has been re- 
moved by an uniform magnetic field B applied in the 
direction of propagation of the applied electric field. The 
levels \e±) are shifted about line-center by an amount of 
=F^B_B/ft (/is=Bohr magneton). The field E is detuned 
from the line center by an amount 5 = ujj rg (B = 0) 



— — 1(1 ± i)e 27rlMx+ + (1 =F i)e 27 " Mx ~ | 2 (8c) L0 +g{B = 0) being the atomic transition frequency in ab- 



|e±4| 2 - ^exp [-47rfcZ.Im( X ±)] , (8d) 



sence of the magnetic field. 

The susceptibilities of the a± components inside the 
medium can be written as 



where, I — \£\ 2 is the input intensity of laser field. 

Note that all the measured quantities defined by 
Eq. JSJl are functions of the frequency of the exciting field. 
If the exciting field is spectrally impure, then the Stokes 
parameters (S a ) are to be obtained by averaging over the 
spectrum S(u>) of the laser field. Thus, Ps in Eq. © are 
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where, 27 = 4|d| 2 w 3 /3fic 3 is the spontaneous decay rate 
of the levels |±), \d\ is the magnitude of the dipole mo- 
ment vector for the transitions |±) <-> \g), N is the atomic 
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FIG. 2: Variation of MOR angle 9 with optical density a 
for s = 27 and different laser line- widths 7 C = O.I7 (solid 
line), 7c = O.57 (dashed line), 7c = 7 (dot-dashed line), and 
7 C = 27 (long-dashed line). We have chosen A = 422.67 
nm corresponding to Ca 40 ISo <-> lPi transitions. Note the 
nonlinear dependence of 9 on a for larger 7 C . 



number density, and s = HbB/K is the Zeeman splitting 
of the excited levels. Using these x±, we can now write 
the field amplitude from Eq. (JSJ as 



So = [i+S+e l 



.£- 



(12) 



where, a = 4TrklN\d\ 2 ftvy = {3X 2 /2tt)NI is the optical 
density of the medium. 

In what follows, we will assume that lui = uj +g (B = 0). 
We calculate (ift)'s using Eqs. (|10fl and © numerically 
for different values of -f c and s. We show the results 
in Fig. [21 We clearly see that for j c -C s, the rotation 
angle 9 is linearly proportional to a. But for "f c > s, this 
variation deviates from linearity in large a domain. This 
behavior can be explained in terms of the off-resonant 
components which dominate for the large j c and large a. 

In order to understand the numerical results, we first 
consider the limit of small optical densities whence 



Si = \[2 + aRe{i(p+ +p_)}] , 
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Clearly, the absorption does not contribute to the ro- 
tation angle. We have again recovered the linear depen- 
dence of 9 on a, provided the approximation (|14l) is valid. 
Thus, any departure in linearity of 9 with respect to a 
indicates the breakdown the approximation i|14fl . The 
numerical results of Fig. clearly show the breakdown of 
the mean field description obtained by replacing x's by 
their average values. 

From the Eq. I|lrj|) . we readily see that in the low a- 
domain, by increasing s (or 7 C ) while keeping 7 C (or s) 
constant, the slope of 9 with a decreases. This is clear 
from the numerical results of Fig. [2J Also, for larger 
values of j e , variation of 9 with a is deviated from lin- 
earity. Linear variation of 9 with a is attributed to the 
monochromatic laser field. If the electric field is spec- 
trally impure, then the off-resonant components also con- 
tribute to 9, through the relations ©, @, and 10. Thus, 
9 starts varying with a linearly in low a limit, then satu- 
rates, and finally decreases to zero to change the direction 
of rotation, for larger 7 C (see Fig. El) • But for smaller j c , 
the linear behavior is retained even for larger a, as the 
off-resonant components are not dominant in this param- 
eter zone. 

Next we consider the variation of the degree of polar- 
ization P and the ellipticity e with a. We have noticed 
that P decreases from unity for increasing a. This means 
that the output field no longer remains fully polarized, 
rather it becomes partially polarized. 

Again, from the Eqs. I|6d(l . (|8d|l . and ((TIT) , it i s clear 
that an integration over the entire range of detuning 5 
would yield (S3} — 0, as the integrand is an odd function 
of oj. Thus the ellipticity e becomes zero. This means 
that the polarized part of the output field remains linear. 

From the above discussion, it is clear that the output 
field gets rotated as a manifestation of cumulative effect 
of optical density, magnetic field, and laser line-width. 
Also it becomes partially polarized with no ellipticity. 



Thus, the departure of the rotation angle from the lin- 
earity has to do with the averages of the exponentials 
appearing in I's [Eq. JHJ]. If one were to make the ap- 
proximation of replacing all x's in Eq. © by their aver- 
ages, i.e., 

(exp [2iriklx±]) = exp[27rifc^(x±)] , (14) 
then the Stokes parameters (5*1) and (6*2) would be 
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where, (p±) — (p^) + ((pf) and thus 



IV. QUANTITATIVE MODELING OF 
EXPERIMENTAL RESULTS OF LABEYRIE et al. 
FOR MOR IN SPECTRALLY IMPURE FIELDS 

We now extend our understanding of resonant MOR 
as described in the previous section to explain the exper- 
imental data of Labeyrie et al. In their experiment, a 
cold atomic cloud of Rb 85 is subjected to a static mag- 
netic field. The laser probe beam passing through the 
medium in the direction of the magnetic field is tuned to 
the D 2 line of the atoms (2S 1/2 <-► 2P 3/2 ; A=780.2 nm). 
They have measured the intensities of outputs with differ- 
ent polarizations, as function of laser detuning and also 
at different values of optical density. They have found 
a nonlinear dependence of the MOR angle 9 on optical 
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FIG. 3: Level diagram for the F e = 4 <-> F g = 3 transition. The numbers at the top of the figure indicate the magnetic 
quantum numbers of the sublevels. The relevant Clebsch-Gordan coefficients for the corresponding transitions are given by 
ai = -1/V42,a 2 = -\/§/Z\/\A,a 3 = -1/2^21, a 4 = - v / 5/2v / 21, as = -1/6-/7, a 6 = -1/2a/3,o 7 = -1/3, a 8 = -1/6, a 9 = 
— l/v21,aio = — \/TE/6\/7, an = 2/3V7. The Zeeman splitting of the various sublevels are not shown. 
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FIG. 4: Variation of magneto-optical rotation angle 9 with a for magnetic field (a) 2 Gauss (= 27r x 2.8 MHz) and (b) 8 Gauss 
(= 2tt x If. 2 MHz) for laser line-widths 2-y c = 2vr x 0.5 MHz (solid line), 2j c = 2tv x 1 MHz (dashed line), 2j c = 2tt x 3 MHz 
(dot-dashed line), and 2-y c = 2n x 5 MHz (long-dashed line). The dot-dashed curves correspond to the width of the laser used 
in the experiment Q. Note that the line- width of the D2 line is 2-k x 5.88 MHz. 



density. They found for larger magnetic field that the 
linear behavior is recovered. 

To explain these observations, we consider the rele- 
vant energy levels of Rb 85 as used in the experiment 
(see Fig. |3J). The x-polarized electric field © is applied 
to cold Rb 85 medium near resonantly. The medium is 
subjected to uniform magnetic field B applied in the z- 
direction, i.e., along the direction of propagation of |JS}. 

A. Calculation of \± an d optical density 

The circular components a± of the input electric field 
P) l interact with the transitions m e <-> m g = m e — 1 and 
m e <-> m g — m e + 1, respectively. We assume that the 
electric field is weak enough so that it is sufficient to use 
the linear response of the system to the laser field. We 
neglect the ground-state coherences. As we are consider- 
ing the cold atoms, we neglect the collisional relaxations 
and Doppler broadening of the sublevels. We also assume 
that the atomic population is equally distributed over all 



the ground sublevels. 

Using all these assumptions, we can write the suscep- 
tibilities x± f° r the a± components as the sum of the 
susceptibilities of all the relevant m e <-» m g transitions 
in the following way: 

^ _ y-" 1 N \dm c ,m g -e±\ 2 i 

7 ft r< m e .m a ~t~ i\0 -\- S m m ) 

(17) 

where, s misimg = (g g m g — g e m e )s is the relative amount 
of Zeeman shift of the excited sublevel m e with respect 
to the Zeeman shifted ground sublevel m g , g g = 1/3 and 
g e = 1/2 are the Lande-g factors of the ground and ex- 
cited levels, respectively. The factor 1/7 comes into the 
expression (fPTjl as we have assumed equal population dis- 
tribution in all the (2F g + 1) = 7 ground sublevels. The 
coherence relaxation rate r mcjrrlg in Eq. I|17|) is given by 

rra e ,m, g ^ y Tfc,m e 1 (1^) 
k 
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FIG. 5: Variation of degree of polarization P with a is shown for magnetic field (a) 2 Gauss and (b) 8 Gauss, for laser line- widths 
27 C = 2tt x 0.5 MHz (solid line), 2-y c = 2tt x 1 MHz (dashed line), 2-y c = 2tt x 3 MHz (dot-dashed line), and 2^ c = 2vr x 5 MHz 
(long-dashed line). The dot-dashed curves correspond to the width of the laser used in the experiment 



where, 7$ j is the spontaneous relaxation rate from the 
sublevel j to i. Here we have assumed that there is no 
spontaneous relaxation from the ground sublevels. The 



terms cJ TOe , m _ and r m m 's can be calculated from the 
relevant Clebsch-Gordan coefficients (see Fig. |3} 6]. The 
Einstein's A coefficient for the D2 line is known to be 
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where, ( || D \\ ) represents the reduced matrix ele- 
ment of the dipole moment vector d metm . The three 
symbols 3/2, 5/2, and 4 correspond to J, /, and F val- 
ues respectively of the upper levels. Thus all T m<i , m 's in 
(|T7jl are found to be equal to (4o, 3 /3/ic 3 )|(§, §, 4 || D \\ 



i|,3)| 2 /2. 



We calculate the optical density a of the medium, when 
the input light field is resonant with m e = «-> m g = 
transition (5 = 0) in the absence of any magnetic field 
(B = 0). For this, we first obtain the total output in- 
tensity from Eq. (|(jaf) averaged by a very narrow laser 
line-shape, i.e., in the limit 7 C — * 0. Using Eq. I|17fl . we 
thus find that the transmittivity of the medium becomes 
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+ o = j- S \s=o = exp(-a) , (20) 



where, a = (3/7).(3A 2 /2tt)7V/. It should be borne in 
mind that it is different from the definition in Sec. IIIII 



the Faraday angle 9 varies with the optical density a for 
different values of j c and B. Clearly, for j c -c s, the 
rotation angle 6 varies linearly with a. But for larger 7 C 
(> s), the variation of 9 with a deviates from linearity in 
large a. This is because the off- resonant components con- 
tribute to the output intensity. Also, note that for a given 
value of 7c, if s is increased, the linearity is maintained 
even in the large a-domain. This is because for larger s, 
the off-resonant components do not contribute much to 
the output intensity. The resonant frequency component 
is always dominant in the optical density range consid- 
ered. We also note that, as j c increases, the linear slope 
of 9 with a decreases in the small a domain. 

In Fig. we show the variation of degree of polariza- 
tion P with a for various values of the B and 7 C . These 
results reveal that with increase in a, the degree of polar- 
ization deviates from unity, i.e., the output electric field 
not only rotates in polarization, but also it becomes par- 
tially polarized. However, the ellipticity of the output 
field still remains zero as we have argued in Sec. IIIII 



B. Discussions 

Using the above expressions of x± [Eq. (fTTfl ] and 
Eq. ©, we calculate the averaged intensities (In) in dif- 
ferent polarization directions. The Stokes parameters S a , 
degree of polarization P, and the Faraday rotation 9 are 
calculated using the relations (7J. In Fig.0] we show how 



V. CONCLUSIONS 

In summary, we have given a quantitative analysis of 
magneto-optical rotation of spectrally impure fields in 
optically thick cold Rb 85 atomic medium. We have shown 
that the dependence of rotation on the optical density of 
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the medium deviates from linearity due to the finite laser results of Labeyrie et al. 
lincwidth. Using our model, we explain the experimental 
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